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In the presence of a short-distance cutoff, the choice of a vacuum state in an inflating, non-de 
Sitter universe is unavoidably ambiguous. The ambiguity is related to the time at which initial 
conditions for the mode functions are specified and to the way the expansion of the universe affects 
those initial conditions. In this paper we study the imprint of these uncertainties on the predic- 
tions of inflation. We parametrize the most general set of possible vacuum initial conditions by 
two phenomenological variables. We find that the generated power spectrum receives oscillatory 
corrections whose amplitude is proportional to the Hubble parameter over the cutoff scale. In order 
to further constrain the phenomenological parameters that characterize the vacuum definition, we 
study gravitational particle production during different cosmological epochs. 



I. INTRODUCTION 

Inflation ij causally explains the origin of the super- 
horizon density perturbations that seed the structures we 
observe in the universe. During inflation, the physical 
size of any given perturbation mode grows faster than 
the Hubble radius. In particular, observationally rele- 
vant modes start well within the Hubble radius, cross it 
and "freeze" till they reenter the Hubble radius and de- 
velop into galaxies, clusters, etc. 0. Causality within 
the Hubble radius naturally allows us to determine their 
initial amplitude by postulating that perturbations start 
in their vacuum state. 

In Minkowski space there exists a well-defined, unique, 
vacuum state, but in any expanding universe — and in 
particular, during inflation — the notion of a vacuum is 
ambiguous 0]. In de Sitter space there exists a con- 
crete set of vacuum states invariant under the symmetry 
group of the spacetime. Such vacua, the Bunch-Davies 
vacuum and the a-vacua, have received widespread at- 
tention in the recent literature 0. However, strictly 
speaking, an inflating spacetime is not de Sitter space- 
time, de Sitter space admits a time-like Killing vector 
field, whereas the non-existence of a time-like Killing vec- 
tor field is what singles out cosmological spacetimes. A 
Friedmann- Robertson- Walker universe whose scale factor 
grows quadratically in cosmic time (a oc t 2 ) is an inflating 
spacetime, and it is far from obvious why de Sitter space- 
time, where a grows exponentially with time (a oc e Ht ) 
shall be a good description of such a spacetime. More 
importantly, in de Sitter space scalar cosmological per- 
turbations are "ill-defined" [f|, as signaled for instance 
by the divergence of the power spectrum in the de Sitter 
limit. Thus, during cosmic inflation one still has to face 
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the problem of defining a vacuum state in an expanding 
universe that does not admit any particular additional 
symmetry. 

Recent advances in the formulation of quantum field 
theory in curved spacetimes @ suggest that such a pre- 
ferred vacuum state might not exist. Whereas one can 
construct mathematically well-defined quantum field- 
theories in globally hyperbolic spacetimes, they do not 
single out any particular quantum state. Allowed phys- 
ical states only have to satisfy the Hadamard condition, 
which is a condition on the ultraviolet (fc — > oo) behavior 
of two-point functions. In this work we assume physics 
to be unknown above cutoff energies A. The presence of 
this ultraviolet cutoff thus renders the Hadamard con- 
dition inapplicable. Therefore, i) in curved spacetimes 
there is no single preferred quantum state and ii) even 
though there is a class of preferred states, in the pres- 
ence of an ultraviolet cutoff A the condition that singles 
out those states is not applicable. 

Nevertheless, it is still possible to single out a set of 
"reasonable" vacuum states by appealing to the flatness 
of an expanding universe at short distances. In the limit 
where the physical length A of a mode is infinitely smaller 
than the Hubble radius H~~ 1 , the expansion is completely 
negligible. Therefore, in the conventional treatments of 
inflation, the vacuum for a particular mode is chosen 
when X/H^ 1 — > 0, i.e. when the vacuum state agrees 
with the Minkowski vacuum. However, if physical laws 
are unknown below certain cutoff length A -1 , the limit 
A — > leads into a region where physics is unknown 0. 
In the presence of a cutoff, in order to minimize the am- 
biguity related to the expansion of the universe, and in 
order to avoid the region of unknown physics, the best 
one can do is to define the vacuum at the time the phys- 
ical length A of the mode equals the fundamental length 
scale of the theory, A -1 y|- In that case, the conven- 
tional predictions of inflation get modified due to the 
finite effects of the expansion of the universe Q. These 
corrections can be expanded as a series in H/A. Let us 
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stress though, that the corrections we are talking about 
have little to do with trans-Planckian physics (see 01 
for a review). Rather, they arise from the uncertainties 
related to the definition of vacuum in an expanding uni- 
verse. If we knew how to uniquely define the vacuum of 
a field in an expanding universe, there would not be any 
uncertainty at all. In fact, in the limit of no expansion, 
H — > 0, all vacuum prescriptions we consider here agree, 
i.e. to zeroth order in H/A all of them yield the same 
"conventional" power spectrum prediction. 

The precise nature of the small H/A corrections is not 
a question of academic interest only. In some inflationary 
models, H/A might be large enough in order for linear 
corrections to leave an observable imprint on the CMB 
spectrum |llL Il2j | , whereas quadratic corrections are ex- 
pected to be unobservable [13(. In addition, the defi- 
nition of a vacuum also directly affects the number of 
particles produced due to the expansion of the universe. 
By imposing observational constraints on the amount of 
the gravitationally produced particles one can thus gain 
information about the realized vacuum 0, Il4j . 

The goal of this paper is not to assess which vacuum 
choice is the correct one, as we have argued that there 
might be no answer for that question. Our goal is rather 
to take a phenomenological approach and find out what is 
a reasonable set of possible vacua, and how and to what 
extent these possibly different vacua alter the conven- 
tional predictions of inflation. Then, instead of relying 
on theoretical arguments to single out the vacuum that 
was realized during inflation, we shall rely on observa- 
tions to put constraints on possibly realized vacua. 



II. FORMALISM 

For completeness, we shortly review in this section the 
formalism of the generation of perturbations in an inflat- 
ing spacetime. The reader might want to skip to the next 
section and eventually refer back for notational details. 
We also summarize our notation in Tabled 

We consider in the following power-law inflation in a 
spatially flat FRW universe, 



ds 2 = a 2 {ri){djf 



dx 2 



(1) 



Power-law inflation is simple enough to allow a straight- 
forward treatment of the equations of motion of the per- 
turbations, and it is general enough to accommodate 
a wide realistic class of inflationary behaviors. During 
power law inflation the scale factor is given by 



a cc \rj\ 



(2) 



where p > 1. In the limit p — ► oo one recovers de Sitter 
space, a oc —l/rj. Note that conformal time is negative 
during inflation (p > 1). Equation J5J also applies to any 
stage of power-law expansion (0 < p < 1), when confor- 
mal time is positive. For p = one recovers Minkowski 
space. 



The behavior of scalar and tensor perturbations during 
inflation driven by a single scalar field can be described in 
terms of a single scalar variable v(r], x) For scalar 

perturbations the variable v is a particular combination 
of metric and scalar field perturbations, whereas in the 
case of tensor perturbations, the variable v is simply pro- 
portional to the amplitude of the gravitational waves. 
The dynamics of v are determined by the quadratic La- 



grangian 



L = - / drx 
2 



1 dv dv a" 
dx 1 dxi a 



(3) 



where a prime denotes a derivative with respect to con- 
formal time and i,j run from 1 to 3. 

The classical variable v can be quantized following the 
standard rules 0- Upon quantization, v turns into an 
operator v, which can be expanded in Fourier modes, 



d 3 k 



(vk^e^ak+vU^e-^al). (4) 



V2J (2^) 3 / 2 

The mode functions v k {rf) obey the differential equation 



a 
a 



v k = 0, 



(5) 



and the operators a k can be interpreted as annihilation 
operators, [ak, dy] = 5(k — k'), if the Vk satisfy the nor- 
malization condition 



v k v k ~v k v k = -2t. 



(6) 



In terms of creation and annihilation operators, the 
Hamiltonian of the system @ is given by 

H = ^Jd 3 k [{v' k 2 + uj 2 v 2 ) a k a_ k + (7) 
+{\v k \ 2 +Jt\v k \ 2 )ala k + h.c], 
where we have defined the squared frequency 

-i = k 2 --. (8) 

a 

Given a mode expansion, i.e. given a particular set 
of annihilation operators a k , the vacuum of the field |0) 
is implicitly defined by the condition a k \0) = 0. Since, 
in the absence of further requirements Eq. Q does not 
suffice alone to uniquely determine the mode functions 
v k , infinitely many definitions of vacuum are possible |3j . 
In the next section we shall formulate additional criteria 
that constrain the possible vacuum state choices. 



1 For scalar perturbations the quoted Lagrangian only applies dur- 
ing a stage of power-law inflation. See Q, Eq. (11.1), for the 
correct expression in an arbitrary background. 
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Once a vacuum state has been determined, it is pos- 
sible to study the imprint of the vacuum upon observ- 
able quantities, such as temperature anisotropics in the 
cosmic microwave background. For scalar perturbations, 
the amplitude of such fluctuations can be characterized 
by the power spectrum of the Bardeen variable 0] 



c = 



(9) 



where Mpi is the Planck mass M Pl = G^ 1 . The power 
spectrum V is then implicitly defined by the correlation 
function of the £ variable, 



dk sin(fcr) 
k kr 



(10) 



Substituting the mode expansion m to the last defini- 
tions one arrives at the following expression for the power 
spectrum, 



P 



k 3 \v k \ 2 



7TM pl 



(11) 



The value of the power spectrum for a given A: is a 
measure of the mean square fluctuations of the variable v 
over comoving distances r ss 1/k [2j. Hence, the physical 
extent of a field fluctuation labeled by k is 



A 



a 
k' 



(12) 



which depends on time for fixed k. In an expanding uni- 
verse there is a natural physical distance scale given by 
the Hubble radius H -1 , where H = a' /a 2 is the Hubble 
parameter. We denote by 9 the (dimensionless) ratio of 
the physical distance associated to a given mode k and 
the Hubble radius. 



A 



P 



1 



H- 



1 — p kr] 



(13) 



In the limit of no expansion, p — > 0, 8 tends to zero. 
The parameter 9 will play a crucial role in our further 
discussions. 



III. VACUUM CHOICE 

Several vacuum choices have been proposed in the lit- 
erature, and we discuss some of them in the Appendix. 
Here we shall pursue however a different approach. In- 
stead of declaring one of the several different vacua to 
be the correct one, we shall abstract the property all the 
different vacua share and use that property to define a 
general class of "sensible" vacuum states. 

Specifying a set of mode functions Vk (f?) is tantamount 
to choosing a vacuum state |0). Because Vk is determined 
by the equation of motion (J5J , all that has to be done to 
define a vacuum is to specify initial conditions for Vk- We 



TABLE I: Summary of notation 



Symbol 


Meaning 


Equation 


f3 and a 


Bogolubov coefficients 




C 


Bardeen variable 


m 


V 


Conformal time 




9 


Mode length over Hubble radius 




A 


High energy cutoff 




A 


Physical size of a mode 


cm 


V 


Index of Hankel function 




a 


Scale factor 




k 


Labels a perturbation mode 




V 


Power spectrum 


est 


y 


Power-law expansion exponent 


© 


V 


The quantization variable 


© 


Vk 


Mode function 


© 


X and Y 


Vacuum parameters 




Subscript 


Initial time 






d/dr) 





shall require of any vacuum state that in the limit where 
the physical length of the perturbation is much smaller 
than the Hubble radius (6 <C 1), the mode functions ap- 
proach the ones in Minkowski space Vk ~ e - lkrl /Vk. In 
fact, in such a limit cosmic expansion should be irrele- 
vant. In other words, when 9 — > the values of Vk and its 
time derivative respectively approach (up to an irrelevant 
phase) 



1 



Vk = —j= and v' k — —iVk. 



(14) 



All the vacuum choices we discuss in the appendix share 
this property. However, in an inflating universe (where 
H remains finite) physics is not well defined in the limit 
6^0. In that limit the physical length of the per- 
turbation is infinitely smaller than the Planck length. 
Hence, one cannot rely on our conventional understand- 
ing of physics in that regime 2 . Thus, for a given mode k, 
the best one can do is pick a vacuum by prescribing the 
values of Vk and v' k at a finite time rjo such that the phys- 
ical length of the mode is much smaller than the Hubble 
radius (in order for the expansion to be "unimportant" ) , 
but much larger than the Planck scale (in order for trans- 
Planckian effects to be negligible). We shall take r/ to 
be the time when the physical length of the mode equals 
a given (fixed) length scale A -1 , 



(15) 



2 In Pre-Big Bang |l (1 and Epkyrotic/Cyclic |l7l scenarios, the 
universe contracts from Minkowski space at past infinity. In that 
case the limit 8 — > is physically well defined. For those models, 
the ambiguities we are considering in this work can be avoided. 
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The scale A is the highest possible scale where we can 
trust our understanding of physics. Conventionally it is 
assumed that A is the Planck scale A = Mpi w 10 19 GeV, 
although the cutoff could be as low as A sw 1 TeV [l8). 
Note that the time of cutoff crossing r/o is fc-dependent. 

According to Eq. (|f 5JI . initial conditions for each mode 
are prescribed at the same physical length A" 1 . Hence, 
the hypersurface where initial conditions are chosen is 
timelike. For comparison, it is going to be useful to con- 
sider an additional hypersurface where initial conditions 
for the modes are specified. We shall choose this hyper- 
surface to be the constant time spatial section where the 
Hubble parameter equals the cutoff A, 



H(vo) = A. 



(16) 



In this case 7]q does not depend on k. Note that in the lat- 
ter case, observationally relevant modes might be trans- 
Planckian at the time initial conditions are fixed 0- 

Because initial conditions for the mode functions are 
fixed at a finite time rather that at ?7o = — oo, one ex- 
pects corrections to Eqs. I|14|l due to the expansion of 
the universe. The only dimensionless quantity one can 
construct from the two scales that appear to be relevant 
in the problem — the physical length of the mode A and 
the Hubble radius H _1 — is their ratio 6. Hence, those 
corrections are expected to be a function of 6*o, the value 
of evaluated at the time initial conditions are imposed. 
Since we assume 9q to be small at that time, and because 
those corrections should vanish in the limit 9q —> 0, we 
can expand them in a power series around 9q = 0. We 
shall concentrate on the lowest order corrections (first 
order), i.e. we impose 0] 



Y-X 



6o + 0(9 2 ) 



(17a) 



.(17b) 



TABLE II: Properties of different vacuum prescriptions 



Vacuum prescription 


X 


Y 


Conventional 





i(l-2p)/(l-p) 


Adiabatic > 1st order 








Hamiltonian diagonalization 








Danielsson 


— i 


i 



ITT1 (see the Appendix for details). If X and Y were p- 
independent, one could use the information about vacua 
in spacetimes other than Minkowski to restrict the val- 
ues of X and Y. For instance, if it turned out that the 
Bunch-Davies vacuum is the only consistent vacuum in 
de Sitter spacetime one could directly compute the 
values of X and Y for de Sitter and apply them to non-de 
Sitter spacetimes. However, at a phenomcnological level 
X and Y could be p-dependent, as happens for instance 
in the "conventional" vacuum prescription. 

Our next task is to implement the "generalized" vac- 
uum prescription i|17[l during a stage of power-law infla- 
tion. In a power-law expanding, not necessarily inflating 
universe the general solution of Eq. J5J is 

v k (rj) = hi 1 / 2 [A k H u (\k V \) + B k H*(\k V \)} , (19) 

where the H v is the Hankel function of the first kind 
if conformal time is negative (p > 1) and of the second 
kind if conformal time is positive (p < 1). The index v 
is given by 



1 



v = — 
2 



p-l 



(20) 



We shall determine the values of the time-independent 
complex coefficients A k and B k by imposing the condi- 
tions H17|) on the exact solution l|19|l . For large values of 
\krj\ the Hankel function H u has the expansion 



Here, Y and X are two complex parameters, and <j>i 
and fa are two arbitrary real phases. The normaliza- 
tion condition 10 implies that the phases are the same, 
<Pi = 4>2- Because the power spectrum (lllf) is invariant 
under v k — > e^v/. we shall drop the phases, 4>i ~ 4>2 = 0. 
In the same way, Eq. © constraints the values of Y, 

Re(Y) = 0. (18) 

The variables X and Y are two phenomenological pa- 
rameters that characterize the choice of vacuum to first 
order in 9q. Our approach could be generalized to higher 
orders by including additional parameters, but for our 
purposes it will suffice to consider the lowest order cor- 
rections. All the sensible vacuum prescriptions we are 
aware of can be cast in the form l|17l) . In fact, our point 
is that any vacuum prescription that can be cast in the 
form (|17fl is a sensible vacuum choice. The specific values 
of X and Y for the particular vacuum prescriptions that 
have been considered in the literature are listed in Table 



H V (M) 




Av 1 - 1 



O(0 2 



(21) 



, v 1 

-ifc?7=F«7r (2 + 4 



which is valid for both p > 1 (upper sign) and p < 1 
(lower sign). Plugging the expansion l|21|) into (|19fl and 
matching to Eqs. (|17|) we find (to first order in 8q) 



A k = 
B k = 
<p = 



2 2 



i-p 



± 



7T ( 1 - 2p 
I-P 



(22a) 
(22b) 
(22c) 



For 6*o = 0, A k = \/tt/2 and B k = 0. These are the values 
one conventionally chooses when computing inflationary 
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spectra 2] . For finite 9q one hence obtains corrections to 
these "conventional" results. 

Let us emphasize again that the corrections in pow- 
ers of #o we consider here are not directly related to 
trans-Planckian effects. Our corrections become large on 
length scales of the Hubble radius. Their origin can be 
ultimately traced back to the ambiguity in defining the 
notion of a particle in an expanding universe, when the 
Compton wavelength has a size comparable to the Hub- 
ble radius Q . Unknown trans-Planckian physics only en- 
ter our discussion by preventing us from taking the limit 
6q — ► 0, making those corrections finite (but still small) 
rather than zero. 



IV. IMPRINT ON THE POWER SPECTRUM 

The power spectrum IjlOjl is defined through the vac- 
uum expectation value of the two-point function of the 
Bardeen variable £. Hence, the choice of vacuum directly 
affects the power spectrum. 

Substituting the mode function 119|) into Eq. 
using the values of the coefficients (|2*2*|l and taking Eq. 
(|18|1 into account one easily finds in the long-wavelength 
limit {krj <C 1) 



V = V C 



f - Re(X) 9 cos 



2p6» ( 7 1 +7r(l-2p) 



Im(X) 9 sin 



1-p 
2p6> ~ 1 +7r(l-2p) 
1-P 



■O{0l) ,(23) 



where Vc is the conventional power spectrum generated 
during power-law inflation, 



V c = 



|i»l 2 /2(p-i; 



/ / 



x p 



H 2 

M 2 
m Pl 



(24) 



In the last formula, denotes the value of the Hubble 
parameter when the particular reference mode crosses the 
cutoff A. The spectral index n$ of the conventional power 
spectrum is directly related to p by the equation 



ns = 1 - 



p-1 



(25) 



Current observations |2Ctl2ll | set the lower limit ns > 0.9, 
which implies p > 21. Notice that in the limit p — > oo, 
the conventional power spectrum diverges, Vc oc p. In 
de Sitter space the theory of scalar cosmological pertur- 
bations is not well-defined [5j. 

The value of 6$ depends on the hypersurface at which 
initial conditions are chosen. If the initial time is chosen 
to be cutoff-crossing, Eq. (|15fl . 8q is given by 



?o = 



X 



-1/p 



(26) 



On the other hand, if the time is chosen to be the same 
for all modes, Eq. (|16fl . 6 is given by 



~~k 



(27) 



In both formulas H„ is again the value of H at the time 
a comoving reference mode crosses the cutoff. 

The term in the square bracket of Eq. I|23|l contains 
the corrections to the standard predictions of inflation 
due to the ambiguity in the choice of vacuum. For arbi- 
trary values of X these corrections are oscillatory, with 
amplitude 9q. First order 9q corrections are absent if and 
only if 



X = 0. 



(28) 



Taking into account the constraint (|I8|I . our generalized 
vacuum choice spans a 3-dimensional space parametrized 
by, say, Im(F), Re(X) and lm(X). Equation J2HJl just 
says that the set of parameters for which there are no 
linear corrections is a two-dimensional plane, i.e. it is of 
zero measure in parameter space. 

If the vacuum is chosen at cutoff crossing, Eq. I|26ll . 
corrections are hence generically linear in iJ»/A |il l22l. 
For large values of p, the frequency of the oscillations 
lu osc [as a function of log(fc/£;*)] is mildly fc-dependent. 
To lowest order in it is given by 



p-1 



(29) 



In the de Sitter limit, p — > oo, there are no oscillations. 
A plot of both the conventional and the corrected power 
spectrum is shown in Figure ^ 

If the vacuum for each mode is chosen at the time when 
the Hubble parameter equals the cutoff, Eq. I|27|) the 
amplitude of the corrections is proportional to (H*/A) p . 
Therefore, unless ff, w A those corrections are highly 
suppressed, since for reasonable values of the spectral 
index, p is large. 

In the Appendix we discuss many of the vacuum pre- 
scriptions that have been proposed in the literature and 
how they translate into specific values of X and Y. The 
results are summarized in Table ^ Inspection of that 
table shows that only for the Danielsson prescription 
corrections to the power spectrum are linear in H/A. 
For the remaining prescriptions, the adiabatic vacuum 
(of order bigger than one) and Hamiltonian diagonaliza- 
tion [23l |26| , the lowest order corrections are at the most 
quadratic 3 . Let us note however that in our opinion, no 
strong theoretical argument singles out any of the differ- 
ent vacuum choices. 



3 This conclusion is in agreement with L'.'l , although the adiabatic 
vacuum we discuss here is not the "adiabatic vacuum" considered 
by those authors. 
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FIG. 1: A plot — in arbitrary units — of the conventional 
(dashed) and the corrected power spectra (solid) vs. 
log(fc/fc»). The values of the parameters are p — 50 (un- 
corrected n s = 0.96), H t /A = 10" 2 and 9 = -4+ 8i. 



V. PARTICLE PRODUCTION 

Given the ambiguities in the definition of a vacuum 
state in an expanding universe, one might take a phe- 
nomenological approach to constrain plausible vacua. In 
general, one expects particles to be produced due to the 
changing gravitational field in an expanding universe. 
The amount of particles produced depends on the way 
the vacuum is defined. Hence, by requiring the rate of 
particle production to be negligible, one can put con- 
straints on the vacuum choice [14J . 

Consider a graviton (or a massless scalar field) prop- 
agating in a fixed expanding universe. The amplitude h 
of any of the two polarization states of the graviton is 
described by the action ©, where v = h/a. Hence, we 
can use our previous results to define a vacuum state for 
gravitational waves. Let us focus on a single comoving 
/c-mode (we drop in the following the subscript k). ft 
is reasonable to expect that the vacuum prescription for 
gravitons and scalars are the same. Hence, we shall as- 
sume that the quantum state of the field is the vacuum 
defined by Eq. (jT7Jl at the time of cutoff crossing, Eq. 
(|15|l . The prescription l|17|) fixes the initial conditions for 
the mode functions i>i n , and hence uniquely determines 
the annihilation operators a- m . We shall call such a vac- 
uum the "in- vacuum", <2i n |in) = 0. At an arbitrary later 
moment of time ?7out, an observer would use the same pre- 
scription (|17fl (with the same values of the parameters X 
and Y) to define a new vacuum |out). He would use a 
constant time hypersurface to define a vacuum state at 
that moment. Because 9 changes with time, such a pre- 
scription yields a different set of modes v ou t and hence, a 
set of different annihilation operators a out . Let us call the 
latter vacuum the out-vacuum, a ou t|out) = 0. In general 
the in- vacuum and out- vacuum are different. In particu- 
lar, the in-vacuum contains particles of the out-vacuum. 
It is easy to show that if the out-modes are expressed in 
terms of the in-modes, 



the number of out-particles iV ou t 
the in-vacuum is 

(in|AU t |in) = 



a out a out 



contained in 



(31) 



This phenomenon is known as gravitational particle pro- 
duction. Using Eq. (|3U[) and the solution l|19|l one can 
easily compute the Bogolubov coefficient f3, 



— — (^out-Sin — B out Aj n ). 



(32) 



Here, A and B are the coefficients (|22|l evaluated at the 
corresponding time: at cutoff crossing for the "in" sub- 
script, and at rj ou t for the "out" subscript. Note that the 
Bogolubov coefficient vanishes if ?7 ou t is taken to be the 
time of cutoff crossing; the mode is indeed in the vacuum 
at cutoff crossing. 

The number of particles in each fc-mode can be used to 
estimate the total energy density of particles produced by 
the gravitational field at any arbitrary moment of time 
77 out . The differential particle density dn is (we restore 
the fc-subindex in the following) 



dn 



J_\jh[ 
47T 2 a 3 



■k 2 dk. 



(33) 



Since we are working to lowest order in we can assume 
the energy of the graviton to be given by E — k/a. Thus, 
the differential energy density of produced particles is 



dp - 



J_\£k 
4tt 2 a' 



-k 6 dk. 



(34) 



The total energy density contributed by the created par- 
ticles is obtained by integrating Eq. (|34(l over a suit- 
able range of modes. At any given moment of time 7y ut, 
the maximum value of k is determined by the cutoff, 
fc max = a(rj out )A. Because our vacuum prescription 117|) 
only makes sense for modes well-inside the Hubble ra- 
dius, we shall take k m i n to be determined by the size of 
the Hubble radius at time rj ou t, fc m in = a(?7out) H out - Sub- 
stituting Eq. Ij32(l into Eq. I|34(l and integrating between 



H out 

16tt 2 



A 



\X\ 2 -I 



A 



(35) 



where I(r) is the integral 



7(r > - 1 { 



(36) 



2a; 2 - 1 /P( 



T^(* 1/P -*)]} 



>dx. 



Win = a v out + P v* 



(30) 



For large values of A/H ou t the oscillatory term in the in- 
tegral averages out. Then, for p > 1/2 the integral / is 
dominated by the upper limit (x = 1), i.e. by the high 
momentum modes, and the integral is of order one. For 
values of p smaller or equal 1/2 the integral is dominated 
by the lower limit (x — H out /A), i.e. by the modes close 
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to the horizon. Because for those modes our vacuum pre- 
scription breaks down, and because horizon size modes 
at that time crossed the cutoff during inflation, we shall 
restrict our attention to p > 1/2 (matter domination or 
cosmic acceleration). 

One can use Eq. i|35|) and the Friedmann equation 
H 2 = 87r ( o cr i t /(3Af|> ; ) to compute the ratio of produced 
particles to the total energy density in a flat universe. As 
an example, we shall consider our recent past, where the 
universe has been mostly matter dominated (p = 2/3). 
Then, I(r) m 3/2 and we find 

Presently, several cosmological probes constrain the con- 
tributions of the different universe constituents to the 
critical energy density j^jj. The uncertainties in those 
contributions are within 10%. Thus, requiring the en- 
ergy density of the produced particles to be less than 
1/10 of the critical density, we find 

If one sets A = Mpi one finds \X\ < 1.1, which is still 
consistent with the Danielsson prescription. For lower 
values of the cutoff, condition l|38l) is even laxer. Sim- 
ilar limits apply during inflation, since the inflaton has 
to remain the dominant energy component during that 
epoch. 



VI. SUMMARY AND CONCLUSIONS 

Due to the ambiguity in the choice of a vacuum state 
during inflation, we have followed a phenomenological ap- 
proach to characterize different possible vacuum choices. 
The vacuum ambiguity has two sources. The first one 
is the time at which the vacuum is defined. We have 
considered two alternatives: cutoff crossing and constant 
H. The second source consists on the way a particu- 
lar vacuum is chosen at the specified time. We have 
parametrized any possible vacuum by a set of two com- 
plex parameters X and Y, subject to the constraint 
(|18|l . These parameters characterize the departures from 
Minkowski vacuum due to the expansion to first order 
in a properly chosen expansion parameter: the ratio of 
physical length of the mode to the Hubble radius. 

To first order, only the parameter X enters the pre- 
dicted power spectra generated during a stage of infla- 
tion. Generically, if initial conditions are chosen at cut- 
off crossing, corrections to the conventional predictions 
of inflation have the form of oscillations with an ampli- 
tude proportional to H/A, Eq. By "generically" 
we mean that the set of X values that yield no linear 
corrections is of zero measure in parameter space. Never- 
theless, for the particular vacuum choices that have been 



traditionally discussed in the literature (adiabatic vac- 
uum and Hamiltonian diagonalization) , corrections are 
at most quadratic. If initial conditions are specified at 
a constant time, the amplitude of the oscillations is pro- 
portional to a large power of H/A, i.e. it is highly sup- 
pressed. 

Since theoretical arguments might not be able to con- 
strain X and Y , we have attempted to set limits on their 
values by studying graviton creation in an expanding uni- 
verse. By requiring the energy density of the produced 
gravitons to be less than ten per cent of the critical en- 
ergy density today we obtain a constraint on the value 
of X, Eq. (|38|l . Even if the cutoff is chosen to be the 
Planck energy, the Danielsson prescription is still allowed 
by observations. For lower values of the cutoff, our ob- 
servational constraint is not significant. 

We conclude that, at the level of our analysis, the is- 
sue about the choice of vacuum is observational, rather 
than theoretical. This does not mean however that in- 
flation is not predictive. All the vacuum choices we have 
discussed yield the same results at zeroth order in our 
expansion parameter. The first order corrections are 
"Planck" -suppressed, i.e. they are proportional to H/A. 
Thus, these corrections are always small, albeit in some 
cases they might be big enough to be detectable. 
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APPENDIX A: VACUUM PRESCRIPTIONS 

We discuss in the following most of the different vac- 
uum choices that have been proposed in the literature. 
They are based on imposing different conditions on the 
vacuum state. Up to what we shall call the "conven- 
tional" choice, they all share the same ambiguity related 
to the time where initial conditions for the mode func- 
tions shall be set. 



1. Adiabatic vacua 

The adiabatic vacuum 0, |25| has been claimed to 
be the best available notion of vacuum in a spacetime 
that has no static regions || or any further additional 
symmetry. The notion of adiabatic vacuum is slightly 
technical, and has been the origin of some confusion in 
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the literature. We hence review here the presentation of 
0. 

The first ingredient of the adiabatic vacuum prescrip- 
tion is the exact, formal "positive frequency" WKB so- 
lution of the equation of motion (JSJ , 



Vk 



: CXp 



W k (fj)dfj 



where Wk is implicitly defined by the equation 



2 _i_(wi_i_wl 

k 2\W k 2 Wl 
and u>k is given by Eq. JSJ, 



2p- 1 



(Al) 



(A2) 



(A3) 



The reader can easily verify that the function IjAlfl sat- 
isfies the normalization condition 

In the limit where the expansion of the universe is 
"slow" , one expects an expansion in the number of deriva- 
tives to be a good approximation to solve for Wk- To 
perform such an expansion one formally replaces 77 by 
77T and expands all quantities in a power series in T . 
To nth adiabatic order, only powers up to T~ n are kept. 
Such an expansion is the second ingredient of the adia- 
batic prescription. It can be easily verified that to lowest 
and first adiabatic order lu^ = k and Wk = ^k = k. Non- 
trivial time-dependent corrections do not show up until 
2nd order. 

If one substitutes the nth order adiabatic approxima- 
tion to Wk into IjAlfl . one obtains the nth order adiabatic 

(n) 

mode ■ Note however that this nth order approxima- 
tion is not an exact solution of the mode equation (JSJ 
anymore. The third step of the adiabatic prescription 
consists of matching the nth order adiabatic approxima- 
tion to the exact mode solutions of Eq. JSJ at an 
arbitrary time 770- In other words, the nth adiabatic ap- 
proximation is used to prescribe the values of Vk and v' k 
at time 770, 



v k (vo) = 4 n) M, VfcM = v'k m {rj ) 



tin). 



(A4) 



Because the WKB solution has been expanded only to 
nth adiabatic order, the matching of the two functions 
has to be carried only to that same order, i.e. only terms 
containing up to n powers of 1/T have to be matched. 
Conditions (|A4|) then uniquely define the mode function 
Vk- By construction, Vk is an exact solution of Eq. JSJ, 
with initial conditions determined by the "distorted" adi- 
abatic modes ul at time 770- 

If one carries over the program described above, one 
obtains expressions for Ak and Bk as a power series in 
l/(T?7o). To nth adiabatic order only terms up to order 
l/(Tryo)" are significant. Hence, if we are interested in 
determining Ak and Bk to first order in 8q oc l/(Trjo), it 
hence suffices to consider the 1st order adiabatic vacuum; 



higher order adiabatic vacua only yield higher order cor- 
rections. Because to first adiabatic order Wk = k, the 
matching condition l|A4(l to that same adiabatic order 
reads, up to an irrelevant common phase, 

Vk(vo) = -/=, '"kiVo) = -iVk. (A5) 

By comparison with Eqs. I|17|) we thus find 

Y = 0, X = 0. (A6) 



2. The conventional choice 

Conventionally, the vacuum is chosen by requiring that 
the mode functions Vk reduce to the Minkowski ones in 
the limit 77 — > — 00, 

Ufc(-oo) = -j=, Ufc(-oo) = -iVk, (A7) 

Using the expansion l|21(l one can directly find the val- 
ues of the coefficients Ak and Bk for which the mode 
function ljl"§|) satisfies Eqs. i|A7|l . 



.4, 



2' 



B k = 0. 



(A8) 



Once these coefficients are known one can then determine 
X and Y by comparing them Eqs. I|22l) . 



X = 0, Y 



l-2p 
1-P ' 



(A9) 



The conventional vacuum is sometimes called the adia- 
batic vacuum, since it agrees with the adiabatic vacuum 
(of any order) in the limit 77 — *• —00. One should bear 
in mind however that the conventional vacuum is only a 
particular adiabatic vacuum in the infinite set of adia- 
batic vacua parametrized by 770- Let us note in addition 
that all the vacuum prescriptions we discuss here also 
yield the adiabatic vacuum in the limit ?/o — » —00. 



3. Hamiltonian diagonalization 

The Hamiltonian diagonalization vacuum is the state 
that diagonalizes the field Hamiltonian J7J) at a given mo- 
ment of time 770 |2J,|3j. Alternatively, it can be shown that 
the Hamiltonian diagonalization vacuum is the state that 
minimizes the energy E — (0|W|0) at 770. It directly fol- 
lows from the Hamiltonian (JJJ and the normalization con- 
dition © that the diagonalization prescription amounts 



4 During the preparation of this work the preprint |2(I appeared, 
which discusses how the choice of different Hamiltonians affects 
the state of minimal energy. 
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to imposing the initial conditions 



v k {vo) = , = , v'kiVo) = -W^kivo), (A10) 

where uj k is the frequency © ■ Note that it is only possi- 
ble to diagonalize the Hamiltonian for a particular subset 
of field modes. 

The squared frequency l)A3|l is quadratic in 7/0 • Hence 
the expansion of the square roots in Eq. I|A10|) in a power 
series in 6q does not contain linear terms (oc 60). There- 
fore, we immediately find 

X = 0, Y = 0. (All) 

In the absence of a cutoff, Hamiltonian diagonalization 
has has been criticized on the grounds of excessive parti- 
cle production Q- Note that the Hamiltonian diagonal- 
ization and 1st order adiabatic vacua are the same. 



4. Danielsson prescription 

According to Danielsson 0, , the vacuum shall be 
chosen by imposing the following condition on the mode 



functions, 

(-)' M = -ik-(vo). (A12) 

Using the normalization condition JBJ, Eq. (|A12|I implies 
the following relations (up to a common irrelevant phase) 



v k (r) ) = -j=, v' k (r, ) = -iVk(l + i6 ). (A13) 

Therefore, by comparing Eqs. (|17|l with Eqs. I|A13|) it 
follows 

X = -i, Y = i. (A14) 

The Danielsson vacuum is a state of minimal uncer- 
tainty (23, but not the only one. It is not the state 
of minimal energy either, in the sense that it minimizes 
the Hamiltonian JJJ (see also [29 )■ Nevertheless, at this 
stage the prescription l|A13|l is as valid as any other. 
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